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ABSTRACT 

During the epoch of reionization, local variations in the ionized fraction (patchiness) imprint 
arcminute-scale temperature anisotropies in the CMB through the kinetic Sunyaev -Zel'dovich (kSZ) 
effect. W e employ an improved version of an analytic model of reionization devised in lFurlanetto et alJ 
(2004bj) to calculate the kSZ anisotropy from patchy reionization. This model uses extended Press- 
Schechter theory to determine the distribution and evolution of H II bubbles and produces qualitatively 
similar reionization histories to those seen in recent numerical simulations. We find that the angular 
power spectrum of the kSZ anisotropies depends strongly on the size distribution of the H II bubbles 
and on the duration of reionization. In addition, we show that upcoming measurements of the kSZ 
effect should be able to distinguish between several popular reionization scenarios. In particular, the 
amplitude of the patchy power spectrum for reionization scenarios in which the IGM is significantly 
ionized by Population III stars (or by mini-quasars/decaying particles) can be larger (or smaller) by 
over a factor of 3 than the amplitude in more traditional reionization histories (with temperature 
anisotropies that range between 0.5 and 3 fiK at I = 5000). We highlight the differences in the 
kSZ signal between many possible reionization morphologies and discuss the constraints that future 
observations of the kSZ will place on this epoch. 

Subject headings: cosmology: theory - intergalactic medium - cosmic microwave background 



1. INTRODUCTION 

Scattering of cosmic microwave background (CMB) 
photons off objects after recombination imprints hot and 
cold spots in the CMB. Measurement of these secondary 
anisotropies will elucidate details of the formation and 
evolution of structure in the universe, including the mor- 
phology of reionization. CMB detectors (such as BIMA, 
CBI and soon ACT and SPT 4 ) are beginning to reach 
small enough scales (< 5 arcminutes) where the primor- 
dial temperature anisotropies no longer dominate the 
secondaries. On these scales, the principal secondary 
anisotropy at most wavelengths is the thermal Sunyaev- 
Zel'dovich (tSZ) effe ct, which comes from scatter ing off 
hot intracluster gas l)Zel'dovich fc Sunvaevl ir969'l. The 
tSZ signal is dominated by nearby clusters and therefore 
is not optimal for studying the high-redshift universe. 
However, the unique frequency dependence of the tSZ 
anisotropies (vanishing at ~ 217 GHz) facilitates their 
removal from the signal. This may allow future CMB 
missions to detect the kinetic Sunyaev-Zel'dovich (kSZ) 
effect — owin g to scattering off of ob jects w ith bulk pecu- 
liar motions ijSunvaev fe Z erdovich il980(l . 

During the epoch of reionization, local variations in 
the ionized fraction (or patchiness) contribute to the 
kSZ signal. The anisotropies produced by this patch- 
iness have been cal culated for various analytic mod- 
els of reionization l)Cruzinov fc HiJ 119981 iKnox et all 

1 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., 
Cambridge, MA 02138; mmcquinn@cfa.harvard.edu 

2 Division of Physics, Mathematics, & Astronomy; California 
Institute of Technology; Pasadena, CA 91125 

3 Jefferson Laboratory of Physics, Harvard University, Cam- 
bridge, MA 02138 

4 For more information, 
www.astro.caltech.edu/'- 



http 



-ti 

' / www .hep.upenn.edu / act / act .html 
//astro. uchicago.edu/spt/ respectively. 



see http://bima.astro.umd.edu 
ip/CBl 



11993 IValageas et al.ll20f)Tl ISantos et aLll20fil . These 
treatments differ considerably in their description of 
this epoch, but most find that the kSZ temperature 
anisotropies produced during reionization have a similar 
amplitude to the kSZ anisotropies produced after reion- 
ization. This suggests that observations of the kSZ effect 
might provide important constraints on the reionization 
epoch. 

Existing probes of this era have provided few defini- 
tive answers about reionization. Observations of Lya ab- 
sorption in the spect ra of high-redshift quasars indicate 
that it ends at z ~ 6 feecker et al|l200lHFan et al.ll2002t 
i White et al.ll2003l iSokasian et alJl2003t iWvithe fc Loebl 
l2004UMesinger fc Haimanll2004t lOh fc Furlanett oll 2005) . 
The main difficulty with these measurements is that 
the Lya optical depth is extrem ely large in a neutral 
medium l|Gunn fc Petersor]ll965[) . making it difficult to 
place strong constraints when the neutral fraction ex- 
ceeds ~ 10~ 3 . On the other hand, measurements of 
the cosmic microwave background by the Wilkinson Mi- 
crowave Anisotropy Probe ( WMAP) suggest a high op- 
tical depth to electron scattering, a pparently requiring 
reionization to beg in at z > 14 (K ogut et alJ 120031: 
iSpergel et all 120031) . Unfortunately, the WMAP data 
provide only an integral constraint on the reionization 
history. 

In thi s paper, we employ the analytic model out- 
lined in IFurlanetto et all l|2004bl hereafter FZH04) to 
calculate the kSZ signal for many different reioniza- 
tion histories. Most analytic models of reionization 
are based on the growth of H II regions around in- 
dividual galaxies or q uasars l)Arons fc Wingertl Il972t 
iBarkana fc Loebl 120011) , T his contrasts with current 
state-of-the-art simulations llSokasian et alJ 120031 120041 
ICiardi et al .112001 IFurlanetto et al.ll2004dl). which find a 
relatively small number of large ionized regions around 
clusters of sources. Moreover, in these simulations, reion- 
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ization proceeds from high to low density regions, im- 
plying that recombinations play only a secondary role 
in determining the morphology of reionization; instead, 
large-scale bias dominates. The model we employ pro- 
duces similar reionization histories to those seen in these 
simulations. 

This paper is organized as follows. We outline the for- 
malism for calculating the kSZ power spectrum in SJ21 I n 
we review and improve upon the reionization model 
first developed in FZH04. This model is able to pro- 
duce a range of possible morphologies for reionization. 
We then calculate the kSZ power spectrum for several 
different reionization scenarios, highlighting the observ- 
able differences between them (|0J, and discuss whether 
future observations of the kSZ effect will be able to con- 
strain the history of reionization 

In our calculations, we assume a cosmology with fi m = 
0.3, n A = 0.7, n b = 0.046, H = lOO/i km s" 1 Mpc" 1 
(with h = 0.7), n — 1, and os = 0-9, consistent with the 
most recent measurements l)Snergel et al.ll2D03^ . 

2. THE KINETIC SZ EFFECT 

Thomson scattering of CMB photons off free electrons 
with a bulk peculiar velocity produces a temperature 
anisotropy along the line of sight n: 



AT 



■(h) 



d?7 e 



-t(tj) 



an P n • v, 



(1) 



where a is the scale factor, r(rf) is the Thomson optical 
depth to the scatterer at conformal time 77, v(n, 77) is 
the peculiar velocity of the scatterer, and the electron 
number density is 

n e (ri, n) = n e (r]) Xi(rj) [1 + 6 b (f], n) + 5 x (rj, n)]. (2) 

Here, n e is the average number density of electrons (both 
within atomic systems and free), Xi is the global ionized 
fraction, and 5b and S x are the overdensities in the bary- 
onic mass and the ionized fraction (we reserve d for the 
matter overdensity). 5 The two-point correlation function 
w(6) = (AT(n)AT(n')) is 



w(6) = [a T n e (r) )} 2 / dr]W( V ) / drj' W^fiitij 



x< g y>(v^ {m _ V )2 + {v _ V) 2 )V)> (3) 

where q(x, 77) = v (1 + Sf, + 5 X ), W(rj) — a~ 2 x? e _T , and 
repeated indexes are summed. In equation (jSJ, correla- 
tions between points separated by the vector (A77, Ax) 
are treated as equivalent to points separated by (0, Aa;). 
While this approximation should be adequate for our 
purposes, we substitute ^/(qq')(x, 77) (qq')(x, 77') for 
(qq')(a;,77) to capture better the 77 dependence. 

Simplifications to equation © are most apparent 
in Fourier space: equation Q involves the integral 
J driW(ri) n • q(k) e jk '™( 7)0 ~' 7 ), which suffers severe can- 
cellation for k along the line of sight. Thus, for modes 
much shorter than the thickness of the window function 
and in the flat sky approximation (fc = l/x, x = 770 — 77), 

5 Note that cq. J2J is not a formally rigorous perturbation ex- 
pansion. However, it suffices for our purposes because the second 
order & x &b term contributes a negligible fraction of the total kSZ 
signal. 



the Fourier space version of equation O ) reduces to 
UKaiseHllflflatTaffe fc KaTnioTikowsk^TflflaT 

Q = (*t nMf J W{ V ) 2 P q± (l/x, r,), (4) 

where (qj_(k) • ql(k')) = 2 (2tt) 3 6 3 D (k - k') P q± (k) and 
qj_(k) is the projection of q(k) perpendicular to k, 



q±(k): 



x[^(|k-k'|) + 4(|k-k'| 



(5) 



Here, /j' = k-k' and v(k) = v(k) k where v(k) is the lin- 
ear theory velocity field [v(k) = i k (fa/k) 8 lln (k, 77), with 
/ = d log D/d log a and D the growth factor]. The linear 
theory velocity suffices because nonlinear contributions 
to the th e velocity of objects ar e negligible in the early 
universe. iCoorav fc Chenl l)2002|) estimated the nonlinear 
velocity contribution to the kSZ power spectrum by uti- 
lizing virial velocity arguments for spherically-symmetric 
halos. They found that this contributes about an order 
of magnitude less signal than the kSZ effect from the 
linear theory velocities. At the high redshifts relevant 
to reionization, the impact of nonlinear velocities will be 
even less important. 

To simplify equation Q, note that the three-point 
terms in P q± are negligible because they involve corre- 
lations with just v(k) at one point, and vj_(k) = by 
definition. This leaves three four-point correlation func- 
tions. We can decompose these into a sum of all possible 
pairs of two-point functions alo ng with the c onnected 
fourth moments l)Ma fc Frvll2002HSantos et alJl2003(l 

P,Ak) = \ E [ ^{(l~v' 2 )Pv V (k')Pa b (\k-k<\) 



a.b=8b,8; 

(l-^)k' 



±-p av (k')P bv (\k-k'\ 



|k-k, 

+ J v /( 1 -^ 2 )( 1 -^" 2 ) cos ^ - ^ 

xP SaSbVV (k-k',-k-k" ) k',k"), (6) 

where Ps,s t w is the connected fourth moment and 
{(/>', (j>"} are the polar angles of the vectors {k',k"}. In 
the halo model and in our reionization model, where 
H II bubbles are assumed spherical and correlations be- 
tween objects are handled similarly to those in the halo 
model, Ps,8.w does not depend on <j>'. Therefore the 
contribution from the connected fourth moments vanish. 
iMa fc Frvl ((2002) verify in simulations that t he contribu- 
tion to P q± from Pg b s b w is indeed negligible. iZahn et alJ 
(2005) compute the kSZ for a nearly identical distribu- 
tion of bubbles sizes as we do, but allow for more com- 
plicated configurations of bubbles. The agreement of our 
kSZ signal with theirs (see Fig. and H4. II below"), sug- 
gests that the connected fourth moments involving 5 X are 
also small and that this is not just a byproduct of our 
simplified model for reionization. 

To further simplify equation ©, we drop the cross 
correlation terms P$ bV and Ps xV . Even if Pg,s, Pw is 
comparable to Ps, v Ps,v, equation JBJ is dominated by k' 
near zero at the scales of interest, k > 1 Mpc -1 (noting 
that Ps lin s lin and Ps,s lin scale roughly as fc' -3 for fc' > 
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0.05 MpcT 1 and that v oc 5n n /k' ). The terms involving 
Pg, v are 0(k'/k) and are therefore suppressed. 

This argument suggests a further simplification for the 
large k modes of interest. Isolating the case a = b = Sb 
in the sum in equation ffijl. dropping terms of 0(k'/k), 
and evaluating the angular integral, we find 

1 



P 



q ± ,a=b=S b {k) — V r 



,Ps b s b {k), 



(7) 



l Pw(k) llHull200(l iMa fc Frvll200l . 



,,2 — f fcfdfc 
u rms J 2-7T 2 

contribution to P 9i is the dominant term 
and is 



where 
This P 5bS , 

for uniform rcionization and is sometimes called the 
Ostrikcr-Vishniac effect (OV). Tra ditionally, this effect 
is cal culated using linea r theor y l|Ostriker fc Vishniad 
ITflSfil IDodelson fc JubadlT993ft , IHul ll200(1ft first calcu- 
lated the nonlinear OV effect, demonstrating that at 
/ ~ 10 4 the signal is significant ly increased by nonlin- 
earities. We use the halo model ijCoorav fc Shetbll2002ft 
to construct Pss - which we then employ to calculate 
the term Pg b g b in equation (JJJ. Unfortunately, calcu- 
lating Pg b g b also requires knowledge of the distribution 
of baryons inside halos: at large overdensities, baryonic 
matter does not trace the dark matter. To model the 
baryon distribution, we employ a filter function F(x) 
such that Pg b s b (k) = F(k)Pgg(k), where Pgg is the non- 
linear matter power spectrum. The function F(k) filters 
large k contributions from the dark matter power spec- 
trum, such that the baryonic matter is not as clumpy as 
t he dark ma tter (see Appendix El for details). 

Iff! Emh uses the same method to calculate the OV. 
Our OV signal agre es well with the analytic calculation of 
iZhantr et aLl(|2004ft . which uses a similar method but also 
incl udes the nonlinear velo city field. IMa fc Frvl (2002) 
and iCoorav fc Shethl l)2002ft employ a different method 
to calculate the OV effect, modeling the gas distribution 
in halos with a /? profile p(r) = A [1 + (r/r c ) 2 ]~ 3 ^/ 2 or 
a similar function. Their calculation results in a slightly 
smaller signal. Perhaps a superior method for calculat- 
ing the OV is through hydrodynamic simulations. Such 
simulations must cover a huge dynamic range, resolving 
nearby halos and large scale structure at high redshifts. 
Simulations tend to agree relatively well w ith the predic- 
tions of these analytic methods at high-/ (|da Silva et al.l 
l20rTillWhite et al. 112001 1 Zahn et al.ll2005ft . but they pre- 
dict less signal at I < 5000. Such deviations are less im- 
portant for studying secondary anisotropics because the 
primordial anisotropies dominate at these scales. (Con- 
tamination of the primordial anisotropies and the result- 
ing bias in cosmological parameter estimation is a differ- 
ent matter.) 

The other two terms contributing to P q± involve in- 
tegrals over P vv Pg^g^ - the term normally included in 
calculations of patchy reionization - and the cross corre- 
lation 2P VV Pg b g :c - ignored in existing calculations, but 
possibly significant if the H II bubbles occur in over- 
dense regions as simulations suggest. We refer to these 
two terms as the "patchy" terms. Since the universe 
reionizes quickly in some of our models (see Figure OJ) , 
the assumption implicit in equation J2J - fc -1 is much 
shorter than the distance over which the window function 
changes appreciably - may break down at some relevant 
/. This is not a concern for the OV effect, where only a 
small portion of the signal comes from the reionization 
epoch. Simple estimates suggest that at scales where 



the kSZ effect dominates over the primordial anisotropies 
(I > 5000), equation (@J should be an adequate approx- 
imation. However, to be safe, we do the calculation in 
configuration space (eq. |3J. Additional details concern- 
ing the configuration space calculation are given in Ap- 
pendix El 

IHul l)2000ft showed that the high- / polarization resulting 
from scattering of the primordial and kinetic quadrupolcs 
off perturbations in the density and ionized fraction have 
negligible amplitude (10~ 3 — 10~ 2 /iK). This results 
from the small magnitude of the quadrupole anisotropies. 
The high-/ polarization signal from reionization is thus 
too small to affect upcoming measurements of B-parity 
polarization from gravitational lensing or from gravita- 
tional waves. Therefore, we only discuss the temperature 
anisotropies here. 

3. H II REGIONS DURING REIONIZATION 

Re c ent numerical simulatio ns (e.g., iSokasian et alJ 
120031 120041: iCiardi et all 12003ft show that reionization 
proceeds "inside-out" from high density clusters of 
sources to voids, at least when the sources resemble star- 
forming galaxies (e.g., Springel & Hernquist 2003; Hern- 
quist & Springel 2003). We therefore associate H II re- 
gions with large-scale overdensities. We assume that a 
galaxy of mass m ga i can ionize a mass C TO gai, where £ is a 
constant that depends on: the efficiency of ionizing pho- 
ton production, the escape fraction, the star formation 
efficiency, and the number of recombinations. Values of 
C S 40 are reasonable for normal star formation, but 
very massive stars can increase the effic iency by an or- 
der of magnitude ijBromm et al.ll2001bft . The criterion 
for a region to be ionized by the galaxies contained in- 
side it is then / co n > where / co n is the fraction of 
mass bound to halos above some m m i n . We will nor- 
mally assume that this minimum mass corresponds to 
a virial temperature of 10 K, at which point hydrogen 
line cooling beco mes efficient. In the extended Press- 
Schechter model llBond "eTalllT99Ti lLacev fc Coldll99^ 
the collapse fraction of halos above the critical mass m m i n 
in a region of mean overdensity S m is 



erfc 



So 



V2Ki„-* 2 (m )Z )] 



(8) 



where tr 2 (m, z) is the variance of density fluctuations on 
the scale m, cr^m = a ( m min> z) and S c w 1.686, the 
critical density for collapse. With this equation for the 
collapse fraction, we can write a condition on the mean 
overdensity within an ionized region of mass m, 

5 m > 5 B (m,z) = i c -^(()[4- ff 2 (m, 2 )f/ 2 , (9) 

where K{Q = erf" 1 (l - C _1 ). 

FZH04 showed how to construct the mass function 
of H II regions fro m 5r in an analogous way to the 
halo mass function ijPress fc Sch echtcr 19741 iBond et all 
1991). The barrier in equation © is well approximated 
by a linear function in a 2 , 5b ~ B(m) — Bq + Bi<r 2 (m), 
where the redshift dependence is implicit. In tha t case, 
the m ass function has an analytic expression ijShethl 
1998): 



dln< 



din? 



Bn 



a(m) 



exp 



B 2 {m) 
2a 2 (m) 



(10) 
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where p is the mean density of the universe. Equa- 
tion (Ej gives the comoving number density of H II 
regions with masses in the range (to, to + dm). This 
result is rederived in Appendix IU1 The crucial difference 
between this formula and the standard Press-Schechter 
mass function occurs because Sr is a (decreasing) func- 
tion of to. The barrier is more difficult to cross as one 
goes to smaller scales, which gives the bubbles a charac- 
teristic size that depends primarily on X{ . The one point 
function Xi(z) for the linear barrier is (Appendix Q 



Q(z) = Xi(z) = -e 



-2 Bn B\ 



erfc 



B - B x a 
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min 
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min 



■ erfc 



Bo + B^l 
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min 



^2a 



2 

min 



(ii) 



This equation agrees extremely well with the Q(z) one 
finds with the full barrier (equation [[J}. 

All the quantities we want to calculate are well-defined 
in this model. In addition, since the typical bubble size 
is usually larger than the scale of nonlinearities, calculat- 
ing the desired patchy correlation functions is simply a 
matter of exploring the properties of a Gaussian random 
field. 

3.1. Bias of the H II Regions 

Unfortunately, there was an inaccuracy in the formula 
FZH04 used for the linear bias (their eq. 22) that caused 
them t o underestimate the bi as. (Note that the expres- 
sion in Sheth & Tormcn 2002 is also incorrect.) To com- 
pute the correct bias, we consider how the bubble num- 
ber density n(m) varies with the underlying (large-scale) 
density. This is straightforward for a linear barrier, be- 
cause the barrier remains linear after a translation by the 
origin. 

Suppose we are in a region of linear-extrapolated mat- 
ter overdensity Sr and variance <tr at a smoothing scale 
R; note that the linear-extrapolated matter overdensity 
at redshift z is Sr D(z) (we adopt this notation for S only 
for this section) . The comoving number density of bub- 
bles n(m\dR, o~r) is simply the usual expression, equation 
1)10(1 . with the replacements 

a(m) 2 ->a(m) 2 - a R 2 , 

B(m, z) — > B(m, z) — 5r, 

Bq — > B + Bi (7r 2 — Sr. 

Linearizing the resulting expression and noting that the 
Eulerian overdense region is larger by a factor of (1 + Sr) 
(|Mo fc Whitdll996|) , we find 



S x (m)=S R D(z) 



B 1 a 2 



1 + 



B(m)/a(m) 2 - 1/Bp 
W) 



Bo 



(12) 



where ^(to) is the overdensity of bubbles with mass to 
[note that S x = J dmn(m) S x (m)]. The last term is only 
significant if R is small or if we are so close to reionization 
that B\a 2 B ~ Bo- If we neglect this term, the linear bias 
is thus 

hi ^ 1 i B(m)/a 2 (m)-l/B 

b x {m) = l-\ ^— r . (13) 



Interestingly, for sufficiently small bubbles b x < 0. This 
means that small bubbles become rare in large overdensi- 
ties and more common in undcrdcnsc regions. Physically, 
this occurs because overdense regions are more advanced 
in the reionization process such that the small bubbles 
have already merged with larger H II regions. During the 
late stages of reionization, only the deepest voids contain 
galaxies isolated enough to create small bubbles. 
Finally, we define the average bias 



Q- 



dmn(m) V(m) b x {m), 



(14) 



which is the bias averaged over the bubble filling fac- 
tor. The average bias tends to be quite large throughout 
reionization, typically between 3 < b x < 10. The (in- 
correct) prescription of FZH04 resulted in moderate bias 
b x ~ 1 and thus underestimated the clustering of bub- 
bles. 

Working only to linear order in the expansion of Sr 
can be problematic because the bias is so large. This is 
particularly true for Q near unity. In this limit, even 
small overdensities can lead to the nonsensical result 
xs ~ Q [1 + b x Sr D(z)] > 1, where xs is the ionized frac- 
tion in a region of overdensity Sr. For large Q, only the 
most underdense parts of the universe are not ionized, so 
xg cannot be well-approximated as a linear function in 
Sr; the linear bias tends to overpredict the clustering of 
bubbles. Conceptually, the problem is the same as those 
described in §3.1 of FZH04. The linear bias works much 
better for small Q. Fortunately, for large Q the bubbles 
become so large that they are essentially uncorrelated 
and the bias is unimportant. 

3.2. Autocorrelation of the Ionized Fraction 

In the original derivation of (xi x'j} and (xi 6' b ), FZH04 
used a Poisson model that allowed the bubbles to over- 
lap. However, our expression for n(m) (eq. I10|) does not 
actually allow this to happen: the bubbles are explicitly 
constructed to be the largest contiguous ionized regions 
in the universe. Here, we construct improved versions of 
(xj x'j} where overlap is forbidden. 

Unfortunately, this poses a new set of problems. To 
compute the relevant integrals analytically, we must as- 
sume that the bubbles are spherical. This assumption 
tends to suppress the ionized fraction at distances around 
the average bubble radius. The reason is simple. Imagine 
that all the bubbles are the same size. Then reionization 
would be similar to packing oranges in a crate: there 
is no way to do so without leaving some gaps, because 
the oranges cannot overlap. In reality, the bubbles will 
not be perfectly spherical; some will be elongated and 
able to fill these gaps. The most reliable treatment is 
far from obvious. If we completely disallow overlap, the 
suppression of power at the characteristic radius of the 
bubbles causes significant (and spurious) ringing in the 
power spectrum. The other option is to allow for some 
amount of bubble overlap. Qualitatively, one could think 
of this as allowing the bubbles to deform and fill the gaps. 
However, this will overestimate the signal at scales near 
the bubble radii. Of course this problem becomes worse 
as Q — > 1 and overlap becomes more common. 

Our solution is to allow bubble overlap only when the 
ionized fraction is small. Outside of this regime, we note 
that the characteristic bubble radius is > 5 Mpc when 
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overlap is significant at Q(z) > 0.5. Figure Q] shows this 
explicitly: we plot the characteristic bubble radius (de- 
fined as the peak of Vdn/d\ogR) as a function of the 
ionized fraction at several fixed redshifts. When the bub- 
bles become large, the dark matter correlation function 
£,5,5 (r) is so small on the relevant scales that we can ignore 
the correlations between bubbles regardless of the bias. 
Therefore, we break our solution into two regimes: in 
the first, the bubbles are small and both the one and two 
bubble terms, P\ and P 2 , are important. 6 In the other, 
the bubbles are large and only the one bubble term Pi 
is significant. We have 



(xiXi)(r) = 



(l-Q)Pi(r) 
Pi(r)+P 2 (r) 



Q 2 Q > V, 
otherwise, 



where 



Pi(r) = / dmn(m)Vo(m,r), 



(15) 



(16) 



P 2 (r) — J dm\n(m\) J d 3 ri / dm 2 n.(m 2 ) 

x /d 3 r 2 [l + £(n-r 2 |mi,m 2 )]. (17) 

Here, Vo(m,r) is the volume within a sphere of mass m 
that can encompass two points separated by a distance 
r. The ansatz [1 — Q]P\{r) + Q 2 for the large-bubble 
regime is meant to approximate the case when we can ig- 
nore correlations between bubbles; note that it obeys the 
necessary limits described in FZH04. However, our pre- 
scription does miss the less important large scale limit, 
(xix'j) — Q 2 — > b 2 {56'). This does not affect our calcula- 
tion because (66'} and (vv 1 ) are so small on these scales 
- scales where the primordial anisotropics dominate any- 
way. 

The integrand over ri and r 2 in equation (|1 T|) depends 
on our scheme for overlap, which we outline below. The 
correlation function £(r|mi , m 2 ) is the excess probability 
to have a bubble of mass m\ a distance r from a bubble 
of mass m 2 . For simplicity, we use_the mean bubble bias 
b x throughout, so £(r|mi,m 2 ) = b 2 £sg, and we replace 
|ri — r2 1 with max(r, P1+P2), where Ri(mi) and ^2(^2) 
are the bubble radii. This allows us to take £ss out of 
the volume integrals. 

We, therefore, are left with four integrals for the two 
bubble term Pi : two over mass and two over volume. The 
last two can be done analytically if we assume spherical 
symmetry for the bubbles and specify some condition for 
overlap. The most obvious choice is to disallow overlap, 
but, as mentioned, this leads to significant ringing. An- 
other possibility, which we adopt here, is to enforce the 
following conditions: (1) m\ cannot ionize r 2 , and ra 2 
cannot ionize r\\ (2) the center of m 2 cannot lie inside 
mi, but any other part of m 2 can touch mi, and vice 
versa. To ensure that this scheme does not drastically 
overestimate (xjX^) near the bubble edges, we tune rj ~ .5 
in equation (|15Jl to optimize agreement with the "exact" 
result — the correlation functions for the same model con- 
structed explicitly from a Gaussian random field in a box 

6 In analogy with the halo model, the one and two bubble terms 
are respectively the contributions to the correlation functions that 
arise when the points separated by a distance r are inside the same 
bubble (Pi) and are inside different bubbles (P2). 




Fig. 1. — Characteristic bubble radius (the peak of Vdn/dlogR) 
as a function of ionized fraction Xi(Q for z = 6 (solid), z = 9 (long- 
dashed), and z = 15 (short-dashed). The H II regions become quite 
large for Q > 0.5, which justifies our approximation of uncorrelated 
bubbles at large Q. All of the reionization models considered in 
this paper are reionized well before z = 6, so the upper two curves 
are most relevant. 



of side-length 100ft. 1 Mpc using the method described 
in lZahn et aT] l)2005|) . 

Figure [21 shows such a comparison for a model with 
£ = 40. Our pseudo-analytic correlation functions are in 
excellent agreement with the exact result. If we had not 
allowed for some overlap or if we had set rj — 1 in our pre- 
scription for (xi x'j), the agreement would not be nearly 
as good. We do not expect perfect agreement because 
the exact correlation functions use a top hat filter in real 
space as opposed to the /c-space filtering implicit in the 
extended Press-Schechter approach, which results in a 
slightly different bubble mass function (see H4. 41 below'). 

3.3. Cross-Correlation Between the Ionized Fraction 
and Density 

Our prescription for (xi S' b ) is similar to (xiX^). We 
begin with the analog of equation (23) in FZH04 

(xiS b )(r) =—Xi + J dmh^nh(mh) J d 3 r h u(r - r h \m h ) 

x J dmn(m) J d 3 r b [l + &/,(r h - r b )], (18) 

where n/, is the dark matter halo number density, u is 
the normalized halo profile (which we can approximate 
as a delta function because halos are much smaller than 
bubbles), and £bh is the excess probability of having a 
halo of mass m^ at r/j given a bubble of mass m at rj. 

To evaluate equation 1)18(1 . we break up the halos into 
those within the same bubble as the point rf, and those 
outside of it. If resides inside the same bubble of 
mass m as then we know from H3.ll that ^bh,in = 
nh{rnh\fn) / Uhijnh) — 1- The conditional mass function 
nh(mh\m) can be co mputed using the exte nded Press- 
Schechter formalism l|Furlanetto et alJl2004|) . Thus the 
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Fig. 2. — Comparison of the analytic correlation functions in J|] 
(thick lines) with those calculated from a Gaussian random field 
in a box of si de-length lOOfe -1 Mpc (thin lines), using the method 
described in IZahn et"aTI i2005l) . We compare these functions at 
three times during reionization for f = 40. The solid curves are 



(Xi £•) - 



and the dashed curves are £ x $t = (xi 8',}. Note 



that the two methods are in good agreement. Although, we do not 
expect perfect agreement because the underlying mass function of 
the bubbles is slightly different for the two methods. 



contribution to equation (|18(l is 

Pin{r)= J dmn(m)V {m,r) J dm h ^- n h (m h \m) 



dm n(m) Vo (m, r) [1 + Si 



(19) 



In the last line, we use the fact that the inner integral is 
simply 1 + (5s, the mean overdensity of the bubble. 

If i"h is outside the bubble_ containing i"b, then we again 
use the trick £bh w bh{rrih) b x £,ss(r) (with the mean bias 
inserted for simplicity). Thus the contribution to the 
integral in equation l|18|) from outside the bubble is 



Pout(r)- 



dmn(m)Vo(m,r) + / dmn(m) 



d 3 r b [b x ( 55 {r eS )], (20) 



where the rf, integration is over all bubbles that ionize 
rf, but not r/j and, for simplicity, we evaluate £ss at sep- 
aration r e fi = max[i2(m), r] rather than at r — iv Note 
that the second term in P out cancels the first term in P m . 
As before, the last term becomes problematic as Q — ► 1; 
here, these difficulties arise primarily from our adoption 
of linear bias, which tends to overestimate the clustering 
in this limit. Fortunately, the term P ou t is unimportant 
at large Q because the effective bubble radius is quite 
large. So external clustering can be ignored and the one 
bubble term dominates. We again break our calculation 
into two parts: 



(xiS b )(r) = 



Pin ~ Pi 
Pin ~t~ Pout 



Q 



q >/, 

otherwise, 



(21) 



where Pi is given by equation l)l(i[l. The solution p n — 
P\ is an ansatz for when correlations with the density 
field inside the bubbles dominate (xi S' b ), which happens 
when the bubbles become large. We tune rj' to get the 




Fig. 3. — Evolution of the global ionized fraction for single reion- 
ization models with £ = 12, 40, and 500 (dashed, solid and dot- 
dashed, respectively). The dotted curve assumes that f evolves 
linearly with redshift from f = 40 at z ~ 11 to f = 300 at the start 
of reionization (z ~ 22). 



best agreement with the "exact" (xi S' b ) ; this turns out 
to be when r/' ~ .5. Figure [21 compares our expression 
for (xi S' b ) with the "exact" result. We see that we are 
in good agreement. Because we smooth the density field 
on the scale of the bubbles and do not include the effects 
of non-sphericity, we should underpredict some of the 
small scale power. This underestimate appears to be 
very minor, affecting smaller scales than are relevant to 
our kSZ calculation. 

4. THE KSZ EFFECT FROM REIONIZATION 
4.1. Single Reionization Episodes 

In the simplest scenario for reionization, a single gener- 
ation of sources, such as population II stars, ionizes the 
universe. At each redshift, the ionized fraction equals 
the collapsed fraction times £ (©■ We calculate the kSZ 
power spectrum for models with £ =12, 40, and 500, 
which produce total optical depths of r = 0.085, 0.117, 
and 0.184, respectively (assuming that helium is singly 
ionized along with hydrogen during reionization and is 
fully ionized by z — 3; c.f. Sokasian et al. 2002). In 
addition, to investigate the effect of a longer reionization 
epoch on the kSZ signal, we calculate the power spectrum 
for a model where £ increases linearly with redshift from 
C = 40 when reionization is complete (z ~ 11) to £ = 300 
at the start of reionization (z ~ 22). This model yields an 
optical depth of t = 0.138. Physically, a monotonically 
increasing £(z) could happen if lower mass halos have 
larger escape fractions. The optical depths produced by 
all four models are withi n 2-a of the WMA P best fit 
value of r = 0.17 ± .04 <|Kogut et al.ll2003^ . Figure 
plots Q(z) for the four single reionization models. 

Figure 01 shows the various contributions to the high-Z 
kSZ power spectrum for the £ = 40 model. The medium- 
width solid line is the total l 2 Ci/2ir from patchy reioniza- 
tion - the sum of the contributions in equation @ that 
depend on Ps x s x (dashed) and Ps T s b (dot-dashed)/ The 
contribution from Ps^s^ is larger than the cross correla- 
tion Ps^St at all scales, but the latter is never negligible. 

7 Note that P$_s_(k) <-» (xx')(r) through a Fourier transform 
and similarly for the cross correlation. 
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Fig. 4. — Contributions to the power spectrum for f = 40. The 
three curves with the smallest amplitude are the patchy contri- 
butions from Ps x s b (dot-dashed), Pg x s x (dashed), and their total 
(solid). At larger amplitudes are the OV signal (dotted), the sig- 
nal from the lensed primordial anisotropics (thin dashed), and the 
total power (thick solid). In addition, we show the total signal 
for a u niformly ionized mediu m (thin solid) and the patchy signal 
for the|G ruzinov & Hu 1 1998) toy model with effective bubble size 
R = 10 Mpc (solid with triangles). Both have the same Q(z) as in 
the ( = 40 model. 



The dotted line is the OV power spectrum, which is larger 
than the patchy signal. The bulk of the OV effect comes 
from scattering off halos at times well after reionization, 
and it is therefore not optimal for studying the reioniza- 
tion epoch. In addition, we show the total high-Z CMB 
power spectrum, comprised of both the lensed primor- 
dial anisotropics and the kSZ contributions (thick solid). 
Patchy reionization produces a larger kSZ signal than a 
uniform reionization scenario [xi(yi, z) = x~i(z)] with the 
same global reionization history; for comparison, we plot 
in Figure 01 the power spectrum for uniform reionization 
w ith the same Xj as ( = 40 (thin solid). 

iGruzinov fc Hul l|1998fl calculated the patchy signal for 
a simple and popular toy model where (x^(x) Xi(n')) — 

e~( x ~ x ' l 2R (xi — xj) + x\ and where the cross correla- 
tion contribution is ignored. The discussion in <jH shows 
that the assumption of a Gaussian correlation function 
with constant standard deviation R (the effective bubble 
size) throughout reionization is unrealistic. In our model, 
(xi x'^ir) is functionally closer to a decaying exponential 
than a Gaussian owing to the distribution of bubble sizes, 
and the shape of (i;ij)(r) evolves considerably during 
reionization as the bubbles grow and merge. Assuming 
a constant size results in a larger peak l 2 Ci/2w with the 
signal more concentrated at the specific scale . To i llus- 
trate this point, we plot the IGruzinov fc Hul l)1998l) re- 
sults for R — 10 Mpc and the same ionization history as 
our £ = 40 model in Figure 0] (solid with triangles) . 

Figure [S] shows the patchy signal for the four single 
reionization models. Interestingly, this power differs by 
less than a factor of two between the curves. The slight 
differences in the amplitude owe mainly to the increase in 
density with redshift (the probability for scattering goes 
as [1 + z] 3 / 2 ). This effect is partially cancelled because 
the universe ionizes over a somewhat shorter redshift in- 
terval as ( increases. Of the four patchy curves in Figure 
the evolving-^ model has the largest amplitude, which 
reflects a general theme in our results: the amplitude of 




Fig. 5. — Makeup of the high-i signal after removal of the tSZ 
anisotropics . We plot the patchy signal (thick) and total signal 
(thin) for £ = 12 (dashed), 40 (solid), and 500 (dot-dashed), as 
well as a model with variable f (dotted). We include the 1-a error 
bars for ACT assuming perfect removal of point sources and the 
tSZ si gnal. We a lso plot the kSZ power from the simulation of 
IZahn et alJ 120051) (circles) and compare this to our analytic model 
for a comparable Q(z) (thin solid). We do not expect perfect agree- 
ment with the simulation because the bubble size distribution is 
slightly different. 




Fig. 6. — The OV contribution to the signal for three of the single 
reionization models: f = 12 (dashed), 40 (solid), and 500 (dot- 
dashed). We also plot the OV signal for £ = 40 in the extreme cases 
where the gas traces the dark matter distribution at all redshifts 
(thin solid) and with the linear theory density field (thin dashed). 



the patchy signal depends strongly on the duration of 
reionization. Finally, note that the OV contribution to 
the total power is larger than that of the patchy terms 
(see Figures El and . However, the relative differences 
between the OV contributions are comparable to the rel- 
ative differences between the patchy contributions. Fur- 
thermore, most of the differences between the OV signals 
stem from these models fully ionizing at different times 
and not from differences in the duration of the reioniza- 
tion epoch. 

At the bottom of Figure we compare the patchy sig- 
nal predicted by our analytic model to the kSZ power 
compu ted from a 100ft -1 Mpc simulation of IZahn et alJ 
(2005) (their model B). That paper uses the same model 
for the patchy epoch, as outlined in fJ3 imposed on an 
SPH simulation by directly applying the excursion set 
formalism to the density field of the box. The advantages 



to this approach are that it can allow for non-spherical 
bubbles and includes fully nonlinear bias for the bub- 
bles. In addition, the simulation should more accurately 
account for nonlinearities in the density field. The disad- 
vantages are that simulations encounter sample variance 
issues at large scales and offer less direct physical insight 
into the results. The patchy signal from the simulation 
agrees well with our analytic prediction for a comparable 
Q(z), further justifying our simplifications in ^andOU 

The ionization histories for the L20 and S5 simulations 
of lSalvaterra et aD l)2005|) are similar to the single reion- 
ization histories discussed in this section. The authors 
find that the functional form and amplitude of the patchy 
signal for both models is similar, despite reionization 
ending earlier in the L20 model. This conclusion mirrors 
what we find here: the patchy signal is mostly affected by 
the duration of reionization and not by when reionization 
happens. Also, the peak amplitude of their patchy signal 
(« 1.6 x 10 -13 for both models) i s comparable to what we 
find. Despite these similarities, ISalvaterra et alJ (2005) 
note that their 20h~ 1 Mpc simulation box may suppress 
the largest bubbles (particularly at the sizes predicted by 
our model). Such a bias would reduce the signal at lower 
I. It may be that simulations must employ larger boxes 
to accurately predict the patchy signal. 

Figure El also plots the total signal (kSZ effect plus 
lensed primordial anisotropics) for the four models (thick 
lines) as well as the 1-a error bars for the Atacama Cos- 
mology Telescope (ACT) in the 225 GHz channel. The 
error bars suggest that ACT will be able to distinguish 
between three of the four single reionization scenarios. 
These errors were obtained from the ACT specifications 
of 1.1 arcminute resolution and 2 /iK sensitivity, assum- 
ing perfect removal of point sources and the tSZ effect. 
The error bars are reliable if the frequency dependence 
of point source co ntaminants can indeed be mod eled ac- 
curately. However. iHuffenbereer fc SeliaU (|2004^) showed 
that there will be significant bias in the measurements of 
the kSZ signal if this is not the case. It is uncertain how 
well the properties of point source contaminants can be 
modeled on these small scales. 

4.2. Extended Reionization 

Recently a number of theoretical models have at- 
tempted to reconcile the CMB and quasar data by pos- 
tulating an early generation of sources with high ion- 
izing efficiency (most often because they contain mas- 
sive, metal free stars) along with a self-regulation mech- 
anism that switches to normal star formation with 
a lower io nizing efficiency (e.g.. iWvithe fc Loebll200l 
| Ceril2003l iHaiman fc Holderll2003 ISokasian et allbooi 
2004) . Such scenarios can cause a plateau in the ionized 
fraction or even "double" reionization, in which ionized 
phases bracket a mostly n eutral period (although th e lat- 
ter possibility is unlikely; iFurlanetto fc Loebll2005f) . We 
refer to such scenarios as "extended reion i zation" and 
treat them as described in IFurlanetto et all l|2004aj) . 

We model extended reionization with two generations 
of sources. At some specified redshift, the first genera- 
tion of sources turns off. For example, if the first gen- 
eration consisted of massive, metal-free stars, the natu- 
ral self-regulation condition halts the formation of these 
stars when the metallicity in collapsed objects passes 
some threshold ijBromm et alJ2001all200llMackev et alJ 



2003) or when hard photons moderate H2 cooling (e.g. 
Yoshida et al. 2003a, 2004). Another possibility is that 
photohcating halts structure formation in biased regions. 
In either case, feedback slows down the growth of H II re- 
gions by reducing the effective £ in ionized bubbles. We 
will simplify these conditions by assuming an instanta- 
neous transition at some redshift z* . In reality, the tran- 
sition must be smooth and extended ijFurlanetto fc Loell 
2005) , but for our purposes we need only force the reion- 
ization history to stall. Thus, in our model, the universe 
develops a patchwork of H II regions at z* that have not 
yet overlapped completely. The first generation imprints 
a set of ionized bubbles, within which most of the sec- 
ond generation sources grow (because both appear in the 
same overdense regions). The bubbles grow only slowly 
until the total number of ionizations from the second 
generation become comparable to that of the first; after 
this point the evolution approaches the normal behavior. 
One important consequence of such a treatment is that 
it "fr eezes" a scale into the bu bbles for a long period of 
time l)Furlanetto et alJl2004a|) . More realistic histories 
may allow the scale to vary smoothly even during the 
stalling phase. 

Consider a region of mass m with a first generation 
of sources described by (Ci, m min,i), which turn off at 
redshift z* , and a second generation of sources with 
(C2, nT-miafi)- We set m muit 2 = m min - the mass at which 
atomic cooling becomes efficient. The total number of 
ionizations is simply the sum of ionizations from the two 
generations. The excursion set barrier 5 x (m) will be the 
solution of 



l = (Ci-C 2 )erfc 
+ C2 erfc 



S c - 5 x (m) 



v/2[o- 2 (77i m in,i, z') - a 2 (m, z')] 

\ $c - S x (m) 1 

[ \/2[<7 2 (m m in,2, z) - er 2 (m, zj] J ' 



'22) 



where z' = z* for z < z* and z' = z otherwise. Here, 
the complementary error functions are the fraction of 
collapsed gas above the mass thresholds at the two red- 
shifts. With the new barrier, the formalism from 
carries over without further modification (note that the 
barrier remains nearly linear). Of course, our prescrip- 
tion ignores the effect of recombinations. 8 However, even 
a relatively small number of second generation ionizing 
sources should be enough to halt recombinations within 
the H II bubbles. 

Figures \7\ and [S] plot the ionization histories and angu- 
lar power spectra for three extended reionization models 
whose parameters and optical depths are listed in Table 
The first generation of sources turn off at z* = 18 
when xi ~ 0.5 in models (I) and (II). However, the mini- 
mum mass is ten times smaller in model (II) . This smaller 
minimum mass could arise from efficient molecular cool- 
ing within halos. Smaller minimum masses and smaller 
(1 both result in a somewhat smaller effective bubble 
radius i? ff. Because the amplitude of the signal scales 
rough ly as the effective bubble size R c g ijGruzinov "FIT! 
1998), the peak signal in model (II) is smaller than in 
model (I). However, at smaller scales, smaller bubbles 

8 This is actually how we avoid d ouble reionization even th ough 
our transition is instantaneous; c.f. IFurlanetto fc Loebl i2005l) . 
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TABLE 1 

Parameters and optical depths for the three extended 
reionization scenarios considered in this paper. the 
quantity m min is the minimum mass at which atomic cooling 
becomes efficient. 



Model Ci 


C2 


™min,2 


z* 


T 


(I) 12 


500 




18 


0.144 


(II) 12 


100 


0.1 x m m j n 


18 


0.151 


(III) 12 


500 




20 


0.108 




z 



Fig. 7. — Evolution the global ionized fraction for the three 
extended reionization models in Table Model (I) is dashed, 
(II) is dot-dashed, and (III) is dotted. We also plot the single 
reionization model with £ = 12 (solid) for comparison. 



cause (II) to have slightly more power than (I) . In model 
(III) , the first generation of sources terminates at a high 
rcdshift, when Xi ~ 0.15. This results in a smaller signal 
than the other extended reionization scenarios but still 
significantly more than the models described in M.ll The 
OV signals for the three double reionization models are 
nearly identical because the redshift of overlap is almost 
the same. As a result, the differences in the kSZ angular 
power spectra stem primarily from their patchy contri- 
butions. 

All three extended reionization models result in sub- 
stantially more signal at the relevant scales than the sin- 
gle ionization case £ = 12 (see Fig. |HJ. This should allow 
future experiments to distinguish between these two sets 

of scenarios relatively eas ily. The reionization models 

considered in lSantos et alJ l)2003|) are most similar to the 
extended scenarios considered in this section. They find 
that the kSZ signal peaks around ?p Ca k ~ 2000 with an 
amplitude l 2 Ci/2ir ~ 2 x 10 -12 . This / pca k is at a similar 
location in our extended models, and the amplitude is a 
bit larger than ours. However, our signal falls off more 
rapidly than theirs for I > Z pca k. 

We note that our simplified treatment of extended 
reionization, with an instantaneous transition and few 
recombinations, may affect some aspects of these results. 
However, the large amplitude of the patchy component 
depends primarily on the duration over which the patch- 
iness persists and depends only weakly on our model's 
roughly constant characteristic scale during this period. 
Allowing the bubble size to evolve throughout the stalling 
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Fig. 8. — Thick curves show the power spectrum from patchy 
reionization and the thin curves show the total power spectrum of 
the CMB once the tSZ anisotropics are removed. The solid curves 
are for £ = 12 and the dashed, dot-dashed, and dotted curves are 
for the extended reionization models (I), (II) and (III), respectively. 
We include the 1-cr ACT error bars. 



phase will tend to flatten the power spectrum but will not 
substantially reduce the overall power. 

4.3. Uniform Reionization 

On the scale of our bubbles, even a neutral universe is 
optically thin to X-rays and gamma rays. If sources emit 
a large fra c tion of their energy in h i gh energy photons 
(fOhl l200lt IVenkatesan et~alTl200l iRicotti fe Ostrikerl 
l2003HMadau et alJl2004|) or if there is a particle that de- 
cays i nto high energy photons at high redshifts l)Sciama| 
Il982t lHansen & Ha iman 2004f) , such photons could ion- 
ize the universe more or less uniformly. 

First, consider a scenario where a decaying particle uni- 
formly ionizes the universe, which then recombines fully 
before stars turn on. In this case, only the OV power will 
be affected by the uniform component. 9 However, the 
contribution to the OV signal from such high redshifts 
is negligible anyway. If, on the other hand, the epoch of 
uniform ionization overlaps with the patchy epoch, the 
net effect will be to decrease the amplitude of the patchy 
signal and increase the average bubble size. Both sce- 
narios suggest that observations of the kSZ effect may 
be able to break the optical depth degeneracy between 
ionizations by a (uniform) decaying particle and discrete 
sources that would imprint a patchy signal. Alterna- 
tively, the ionizing sources could emit both ultraviolet 
and high energy photons so that the universe has both 
uniform and patchy components to the ionization field. 
This could occur if quasars or mini-quasars account for 
a significant fractio n of the ionizing photons (but see 
Diiks tra et al.1120041 for limits on such a scenario). 

To model this possibility, we suppose that the IGM 
has a uniform ionization fraction x u (z). On top of this 
lie spatial variations from isolated H II regions. In this 
case, the condition / co n > C _1 (l — x u ) replaces the bar- 

9 Of course, the recombination process will lead to patchiness — 
underdense regions will recombine later. We expect this patchiness 
to be minimal because density fluctuations are so small at the 
redshifts we are interested such that a uniformly ionized IGM is a 
reasonable approximation. In principle, given a decaying particle, 
calculating its effect on the kSZ signal would be a straightforward 
exercise. 
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Fig. 9. — Same as Figure^] but with curves for £ = 40 (solid) 
and for three uniform reionization scenarios with ( = 40 (1 — /i) and 
x u = 40/x/ co n. We show /i = 0.25 (dot-dashed), fi = 0.5 (dashed) 
and fx = 1.0 (thin solid). These models have the same Q(z) as the 
canonical £ = 40 model and therefore the same OV signal. We 
include the 1-cr ACT error bars. 



rier in our model: each galaxy can produce a larger 
ionized bubble with the same number of ionizing pho- 
tons. However, rather than varying from zero to unity, 
we have < xh < (1 — x u ). This damps the fluctuations 
from the bubbles, requiring the rescaling (xix'^} — xf — > 
(1 - x u ) 2 {{xix'j) - xf) and (xi5' b ) -> (1 - x u )(x l S' b ) (this 
increases the importance of the cross-correlation term). 
Otherwise our formalism is unchanged. 

In Figure El we plot reionization scenarios where the 
uniform component is proportional to the total ionized 
fraction, x u — fi £ / co ii [it follows that the patchy fraction 
is x p = (1 — li) C/coii]- We plot results for ( — 40 with 
proportionality factors \± = 0.0 (solid), \i = 0.25 (dot- 
dashed), fi — 0.5 (dashed) and fi = 1.0 (thin solid). The 
rough interpretation is that fx parameterizes the fraction 
of ionizations caused by high energy photons emitted by 
discrete sources. The choice fi = 0.0 corresponds to nor- 
mal patchy reionization and fi = 1 to purely uniform 
reionization. The ionized fraction Xi(z) is the same for 
all of these curves (the solid line in Figure |3J). Note that 
the power spectrum for \i = 0.25 and fj, = 0.5 is sup- 
pressed by a factor of ~ 1.5 and ~ 3.0 from the patchy 
model fi = 0.0. The 1-cr error bars in Figure sug- 
gest that ACT is capable of distinguishing between these 
models. 

4.4. Reionization and Uncertainties in the Underlying 
Cosmology 

Even though the least well-known parameter in our 
model is C, there is enough uncertainty in the cosmolog- 
ical parameters and the mass function to strongly affect 
the morphology of reionization. Beyond the usual uncer- 
tainty in f2 m , n, etc., there is also the shape and ampli- 
tude of the mass function at high redshifts; each choice 
will yield a different reionization history. Simulations 
are our best avenue for ascertaining the correct mass 
function. However, to date most high-resolution simu- 
lations have only investigated epochs well after reion- 
ization. We found two studies of the mass function 
at z > 10. Uang-Condell fc HernauistJ l)2001|) showed 
that the Press-Schechter mass function agrees reason- 




z 



Fig. 10. — Comparison of reionization histories with f = 12 be- 
tween cosmologies that have different mass functions and cosmolog- 
ical parameters. The thick solid is the standard cosmology adopted 
in this paper with the Press-Schechter mass function. Other curves 
show the same model except for erg = 0.8 (thin dot-dashed) and for 
tilt n = 0.95 (thin dashed). We also plot the Q(z) for the Sheth- 
Tormen mass function (thick dashed) and Press-Schechter with a 
top hat real space filter calculated from random walks (thick dot- 
dashed). 



ably well with the mass function from their simulation 
at z = 10. However, their l/i -1 comoving Mpc cube 
only allo wed them to prob e the low-mass end. The 
other was iReed et alJ (120031) . which concludes that the 
iSheth fc Tormenl l|2002f) mass function overpredicts sub- 
stantially the mass function in simulations for z > 10. 
In Figure 1101 we plot the expected reionization history 
Q(z) for several mass functions: Sheth-Tormen (thick 
dashed), the usual Press-Schechter function (thick solid), 
and Press-Schechter with a top hat real space filter com- 
puted from random walks (thick dot-dashed). Surpris- 
ingly, the top hat real space filter gives a significantly 
smaller collapse fraction than the sharp /c-space filter 
implicit in the extended Press-Schechter formalism (by 
about a factor of 0.5). This is because the variance in 
the density field within the top hat window function is 
smaller than the variance within the sharp /c-space win- 
dow function. Since the barrier S c — 1.686 is derived 
from top hat considerations, it could well be that the 
top hat in real space gives a more realistic mass function 
at these high redshifts. 

Several observations of large scale structure predict a 
lower value of ctr tha n the a% — 0.9 from CMB data 
l|Tegmarket all 12004]) . We plot Q{z) for ct 8 = 0.8 
(but with other parameters fixed) in Figure 1101 (thin 
dot-dashed). This value for erg causes reionization to 
end later. In addition, CMB observations are consis- 
tent with the tilt n in the primordial power spectrum 
being slightly less than the Harrison-Zel'dovic h choice 
n = 1 HSnergel et a,Ul2003T:lTegma,rk et al.ll20f)l . Even a 
slightly smaller tilt, n — 0.95 (thin dashed), causes the 
universe to reionize significantly later. 

Neither these different cosmological parameters nor a 
different mass function significantly alters the duration 
of reionization. Because the kSZ amplitude is most af- 
fected by the duration of reionization, all of the reioniza- 
tion histories in Figure ITHI should produce similar patchy 
signals. Ho wever, the amplitu de of the tSZ effect scales 
as — o-g -9 ijSeliak et alJl2001^ and the OV signal scales 
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as — <7g~ 6 (jZhang et al.ll2004|) . Thus, if er 8 is smaller 
than 0.9 these anisotropies will be suppressed, and the 
patchy signal will be more prominent. 

If the dark matter is warm, this will have an even 
more drastic effect on the i onized fraction (Barkan a et al.l 
I200H iKitano fc Lowll2005[) . Warm dark matter will cre- 
ate a cutoff in the transfer function at small scales such 
that the collapsed fraction is significantly smaller. While 
warm dark m atter cannot make u niverses with a high 
optical depth (jYoshida et al . 2003b), this can be supple- 
mented by having a decaying particle ionize the universe 
at higher redshifts ({Hansen fc H aiman 20 04^ . 

Finally, we note that because the OV effect is domi- 
nated by lower redshift structure formation, the uncer- 
tainty in the high-z mass function has only a negligible 
effect on that signal. 

5. DISCUSSION 

We find that there is a clear and measurable differ- 
ence between the amplitude of the kSZ signal for models 
where Population II stars dominate reionization, where 
Population III stars are also important (or reionization is 
extended for some other reason), and where there are a 
large number of hard photons (or uniform reionization). 
At scales where the kSZ effect dominates over other pri- 
mordial anisotropies, the amplitude of the patchy signal 
differs by as much as 10/iK 2 . Other aspects of the patchy 
signal also differ between reionization histories, such as 
the multipole I where the patchy power spectrum peaks. 
The peak of the power spectrum is primarily determined 
by the distribution of bubble sizes when Q ~ 0.5 or, for 
our extended reionization scenarios, by the bubble size 
when the ionized fraction "freezes." Ideally, one could fit 
a simple (but physically motivated) two or three param- 
eter model to the observations. Of course, reionization is 
a very complex process and many parameters are needed 
to describe its morphology in detail (in our model, these 
are essentially the parameters determining £ as a function 
of redshift). To add to the difficulty, observations cannot 
separate the OV and patchy contributions to the CMB 
power spectrum, requiring simultaneous fits to both. 10 
The OV power is dominated by post-reionization sources 
and contains much less information about the morphol- 
ogy of reionization than the patchy signal. 

The OV contamination of the patchy signal might not 
be a problem. First, because a substantial fraction of 
the OV power comes from low redshift objects, much of 
it may be removable from the signal. Also, future 21cm 
observations will help separate these two contributions 
(e.g. FZH04. iFurlanetto et al.ll2004at ISalvaterra et all 
20051 iZaldarriaga et all2004D . Furthermore, with a fixed 
cosmology, the OV contribution to the power spectrum 
is determined only by the distribution of gas. To zeroth 
order, we know that gas traces the dark matter distribu- 
tion, which can be successfully modeled with halo theory. 
At large overdensities our understanding of the gas dis- 
tribution is incomplete. But for I < 10 4 , the OV power is 
relatively independent of the gas distribution under rea- 
sonable assumptions. This is evident in Figure El where 
the thin black curve is the OV signal in the most ex- 



treme case in whic h the baryons trace the dark matter. 
iZhang et alJ ((2004) contend that one can reasonably cal- 
culate the OV contribution without any free parameters 
as long as the evolution of the ionized fraction is known. 
The ionized fraction at each redshift is a byproduct of the 
model for patchy reionization. Therefore, the OV signal 
- which has a much flatter power spectrum - is predicted 
by our model with no additional parameters (to lowest 
order) and could help break the degeneracy between the 
redshift of overlap and other parameters. Of course, mi- 
nor differences do remain in the calculation of the OV 
signal between different computational approaches (see 
We expect that future research, particularly sim- 
ulations incorporating realistic gas heating, will help to 
resolve these small discrepancies. 

To this point, we have concentrated on the information 
contained in the kSZ signal. However, it w ill also con- 
tamin a te fits to the primor dial anisotropies . |Kn qx et alJ 
(|1998|) . lSantos et alJ l|2003ft and lZahn et all l|2005fl inves- 
tigated how the signal from patchy reionization will bias 
cosmological parameter determination for WMAP and 
Planck. While the net effect is insignificant for WMAP, 
it will be quite substantial for Planck, which can mea- 
sure C\ to I — 2500. In this case, the bias for many 
cosmologica l parameters i s comp arable to the l-a statis- 
tical error. iSantos et alJ l)2003fl propose adding an ad- 
ditional parameter, the effective amplitude of the kSZ 
power spectrum, in future fits to Planck data, which is 
justified by the fact that the kSZ signal is fairly flat in the 
regime where it affects the fit. A similar conclusion holds 
for our model for reionization (iZahn et alJ 12005'). The 
polarization-polarization CMB power spectrum might be 
more reliable for measuring the cosmological parameters 
at I > 1000 since it is nearly unaffected by patchiness 
l|Zahn et alJl2005^ . 

Other probes of the reionization epoch may help to re- 
strict the set of possible reionization histories, facilitating 
the interpretation of the kSZ signal measured by ACT or 
SPT. A more precise measurement of the optical depth r 
from measurements of large-scale CMB polarization will 
significantly reduce the set of viable models. In addition, 
if future observations of high- redshift quasars or 21 cm 
emission from neutral hydrogen can pin down the redshift 
of overlap, this will again reduce the set of plausible mod- 
els. To the extent that the current large optical depth 
measurement implies a long period of extended reioniza- 
tion, the patchy contribution should extend over a long 
redshift interval and have a substantial amplitude. This 
will allow future experiments to distinguish reionization 
by stellar sources, hard photons, and decaying particles, 
and alleviate many of the degeneracies from large-scale 
CMB polarization data alone. 
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Another "contaminant" of the power spectrum at these scales of the kgz gi j n2 C /2n „ 10 -14 _ 10 -13) th h not obviously 

is the moving clusters of galaxies effect I Molnar & Birkmshawl negligible 
2000 ) . Its amplitude on the scales of interest is smaller than that 

APPENDIX 

THE BARYON DISTRIBUTION 

We construct Pg b g b from the halo model by using NFW dark matter halo profiles with the Jenkins mass function 
and by setting Ps b s b (k) = F(k/kp) 2 Pgg(k), where F(x) filters large k modes in o rder to approximate the e ffects of 
finite gas pressure. The Jeans length is often not the appropriate filtering scale /cp. iGnedin fc Huil l)1998() and lGnedinl 



(Al) 



(1998) suggest kp(z) = 34 y/£l m (z)h Mpc 1 for thermal history-dependent filtering and 

1 



(l + 4x 2 ) lfi 



They find that for S < 10, this filtering function produces errors smaller than 10% in the gas power spectrum Pg b g b . 
We use this scheme for z < 8 and set F(x) = 1 for z > 8, where the gas distribution is less well-understood and it is 
thought that cooling is more efficient. Note that the difference in the kSZ signal between reasonable filtering schemes 
for high redshifts is minute on the scales of interest. 

CONFIGURATION SPACE CALCULATION OF THE KSZ EFFECT 

For the patchy terms, namely (xix'j)(vv') and (x, S' b )(v v'), we perform the kSZ calculation in configuration space 
using equation ©. We first use linear theory to construct the velocity field. In that case, the velocity correlation 
function is 

(v( V , n) vtf, &')) = [£l (r)] n • n' + K,|(r) - U(r)\ (l " &) £ ' &>) , (Bl) 

where v(n) — v(n) • fi, r = x — x' and £|| and £j_ ar e the autocorrelation functions of the velocity parallel and 
perpendicular to r: 

Here primes denote derivatives with respect to kx. If we put the correlation functions constructed here and in Sj3]into 
equation @, we can calculate w(9) directly. For I 1, the angular power spectrum is then 



-f 2 H 2 ol ^kPnnik) (B2) 



Ci » 2tt / 9d9J (W)w(6). (B3) 







DERIVATION OF Q(Z) AND MASS FUNCTION FOR A LINEAR BARRIER 

Our aim in this section is to de rive (or rederive) some relations we w ill need for a linear barrier using the ex- 
tended Press-Schechter formalism ijBond et al.l 119911 ILacev fc Colelll993fl . The advantage of this approach is that 
the probability distribution of 5k - the real-space density within a fc-space top hat filter of radius k - becomes the 
solution to a diffusion equation. This diffusion equation can then be solved for various boundary conditions (corre- 
sponding to the specified barri er) . We first rederive the familiar result for a constant barrier adopting the notation of 
IScannapieco fc Barkanal l|2002|) and then expand the approach to the case of a linear barrier. The variance today in a 
region defined by a sharp fc-space window function is 



S k = a 2 (k) = -L f dkk 2 P lin 
2^ Jo 



(k), (CI) 



wher e Pi n (k) is the linear power spectrum at z — 0. As we change the smoothing scale - or increase k in equation 
(jClll - each step ASk is uncorrelated with the previous Sk (this is not true for other filter choices). Therefore, we can 
write a simple evo lution equation for Q(5, Sk), the probability distribution for the density field at smoothing scale k 
l|Bond et alJll99H) 

Q(6, S k ) = 1 ^ dAS exp f ' A " 

V27T /\bk J-oc 



Q(S-AS,S k ~ASk). (C2) 



2 AS; 

If we expand Q to second order in AS and perform the AS integral, equation (|C2|I reduces to 

Q(S, S k ) - Q(S, S k - AS k ) + l d2Qi6 > S d h p ASk) AS k , (C3) 

noting that (AS) = and that ((AS) 2 ) = AS k . Rearranging the above equation and keeping terms to linear order in 
AS k gives a diffusion equation 

dQ(S,S k ) l d 2 Q(S,S k ) 

dS k 2 OS 2 ' 1 ' 
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In the absence of any barrier, the solution to this equation is 

1 



Qo(S, Sk) 



exp 



2S k 



(C5) 



The solution to a diffusion equation is uniquely specified by the initial and boundary conditions. Thus, for a constant 
barrier B{k) = B$, we can guess the unique solution that satisfies the boundary condition Q(Bq,S) — 0, namely 
Q = Qo($, Sk) — Qn(2B n — 8, Sk)- The second term is an "image " that is placed to cancel the contribution of the source 
term on the boundary ijChandrasekharlfl 943t iBond et al.lll99ll) . 

For the bubble problem, we are ultimately interested in a solution to this diffusion equation for a barrier t hat i s 
linear in the variance B{k) = Bq + BiS k - This amounts to the boundary condition Q(B(k), Sk) = for the PDE (|C4|I . 
To find a solution, let us rescale the variable 5 to y = Bi(6 — BiSk) such that the diffusion equation transforms to 



dQ 

dSk 



Bl 
2 



d 2 Q | 2 dQ 
dy 2 dy 



(C6) 



and the boundary condition becomes Q(B\ Bq, Sk) = 0. Unlike the constant barrier problem, the image method is 
not a frui tful way to solve this problem. Instead, let us look for solutions of the form Q(y,S k ) = g(y) f(Sk)- Then 
equation (|C6|) becomes 



10/ 

fdS k 



2ff 



dg 



X. 



(C7) 



d g 

dy 2 ' dy 

Dropping constant factors, the general solutions for / and g are f(S k ) — exp(XSk) and g{y) = exp[(— 1 ± iX)y] 
where X = — iyl + 2A/B 2 . If we drop the eigenfunctions which do not satisfy the boundary conditions, our solution 
Q(y, Sk) becomes 



!\b 



(y,Sk) 



dXh(X) sm[X(y- B Bx)] exp 



B 



To determine h(X), we impose the initial condition. In 
dition is S D (S) ee limg^p exp[-<5 2 /(2 S k )]/V2^Sj. If 
lims^o exp[-B 2 S k /2-y - y 2 /(2BfS k )}/V2TS k ~ = B ± 8 D {y) exp(- 
Bi/n sin(XBoBi), and we can now integrate equation ()C8|) : 



-y- 

the 
we 

v)- 



i(l+X 2 )5 fc 



(C8) 



(5, S) coordinates, the initial con- 
substitute y for (5 this becomes 
Using this, we find that h(X) = 



Kb 



(y,Sk)- 



1 



V^Sk 
i 

y/2nSk 



exp[ 



B 2 Sk 



y] 





[ y 2 ' 




jexp 


2B 2 Sk_ 


— exp 





r s 2 - 




jexp 


2 Sk_ 


— exp 



|2(B + WBoBxSu) 



(y~2B B^ 
2B\Sk 
8\ 2 



2S k 



(C9) 



(CIO) 



The solution 1C9I) vanishes at the boundary y = BqBi, as required. In addition, equation IjClOl) shows that there is 
an analogous form to the familiar "image" decomposition for the constant barrier problem. 
To obtain the ionized fraction, we integrate over all possible S at S k = <if nin , which yields equation (jlip. We can 



obtain the mass function as well. The probability that a trajectory crosses the barrier in the interval S k to S k 
is 



P(Sk) 



d f BoBl 
dSk J-oo 



dy_ 

Bi 



Q 



lb 



Si dQ 



2 dy 



BqB\ 



Bo 



V^si 



exp 



B{kf 
2 Sk 



(Cll) 



where the second inequality us es the diffusi on equation and the boundary condition Qib{BoBi, Sk) = 0. Equation 
<|C11|) agrees with the results of lShethl lj!998() . which were derived through a different method. The number density of 
bubbles of mass m is then n(m) = po/m P(Sk) \dSk/dm\, which yields equation (JTJJJ. 
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